Dynamical Diagram and Scaling in Polymer Driven Translocation 
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By analyzing the real space nonequilibrium dynamics of polymers, we elucidate the physics of 
driven translocation and propose its dynamical scaling scenario analogous to that in the surface 
growth phenomena. We provide a detailed account of the previously proposed tension-propagation 
formulation and extend it to cover the broader parameter space relevant to real experiments; In 
addition to a near-equilibrium regime, we identify three distinct nonequilibrium regimes reflecting 
the steady-state property of a dragged polymer with the flnite extensibility. Finite size effects are 
also pointed out. These elements are shown to be crucial for the appropriate comparison with 
experiments and simulations. 

PACS numbers; 36.20.Ey, 87.15.H-, 83.50.-v 



I. INTRODUCTION 

Macromolecules can pass through a nanoscopic hole 
which is much smaller in space dimension than the bulk 
coil size. Given the chain-like architecture, this fact itself 
may not be surprising, but the way how it goes threading 
is interesting. Today, this mode of molecular transport, 
called translocation, finds its great potentiality in the 
multidisciplinary field of nano-technology and biological 
sciences e.g., in inventing nano-pore based rapid 

DNA sequencing devices. 

In typical situations, a polymer (with A^o segments) 
passes through driven by the force / acting at the pore 
site. The corresponding setup has been realized in most 
of experiments, where the charged polymers (DNAs) are 
driven by the electric field, and the voltage drop is con- 
centrated in the pore [iHll- From the technological de- 
sign perspective, one of the most important measure is 
the translocation time r(iVo,/), i.e., an average passage 
time for a polymer to go through. For this reason, there 
have been considerable experimental attempts to quan- 
tify the translocation time. [H-O- Extensive numerical 
simulations have also revealed a rich variety of the dy- 
namics Under certain conditions, i.e., when the 
specific pore-polymer interactions can be neglected, these 
results could be described in the light of the universality 
inherent in the chain- like molecules [T^-flsj. It is not, 
however, clear how and to what extent this is possible. 
Diverse "apparent" scaling laws of the translocation time 
reported in literatures remain to be explained (3l-[l6j. 

In this paper, we find three distinct nonequilibrium 
regimes with the different exponents throtigh the tension- 
propagation mechanism proposed in [TtI . fisj . The char- 
acteristic nonequilibrium conformations (trumpet [TP], 
stem- flower [SF], strong-stretching [SS]) of the polymeric 
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FIG. 1: (Color Online) Illustration of various conforma- 
tions during the driven translocation process starting from 
(a) initial equilibrium conformation, (b) Trumpet [TP], (c) 
Stem-flower [SF]. (d) Strong-stretching [SS]. More complete 
processes in respective regimes are shown in Figs. 5-7. The 
tensed "moving domain" is shaded, while the translocated 
part is represented by a dotted curve. 



chain with the finite extensibility are shown in Fig. 1, 
through which we establish a dynamical scaling scenario 
covering various models/situations and broad parameter 
regions relevant to real experiments. 



II. TRANSIENT RESPONSE 

We begin by clarifying the exponents relevant to equi- 
librium (static) and near-equilibrium (dynamic) situa- 
tions. First, the Flory exponent u characterizes the frac- 
tality of the equilibrium polymer coil whose size is given 
by i?o — (a: segment size). Equally important is 

the dynamical exponent z characterizing the dynamics 
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FIG. 2: Dynamical diagram of the driven translocation 
process in the plane of A'o and / (double logarithmic 
scale). Shaded region represents the case to which the near- 
equilibrium theory applies 0, fl^ . fl^ . The arrows (i)-(iii) 
indicate typical sweeping paths of parameters in experiments 
(see Figs. E] and i)). The boundaries /j ~ iVj^", /* ~ 1, 
/** ~ A^o separate regimes characterized by different confor- 
mations (Fig. 1). 



space with excluded- volume/hydrodynamic interactions. 
A drastic simplification is achieved, however, by intro- 
ducing the " two phase picture", i.e., a moving and a 
quiescent domains separated by the tension front, and 
applying the steady-state approximation for the former 
(Fig. 1). It amounts to say that the relevant properties 
of the moving domain can be captured by the represen- 
tative velocity V{t) and the steady-state conformation of 
the N{t) — M{t) segments dragged at its velocity. Leav- 
ing its justification in the Appendix I, we remark here 
that there remains only single gross variable, i.e., the po- 
sition of the front, or the boundary between the moving 
and the quiescent domains. This coarse-graining makes 
the problem analytically treatable. 

Then, by considering the mass conservation relation, 
we can show that the boundary position obeys the fol- 
lowing equation; 
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(1) 



of "critical fluctuation" (polymer coil) To']. Normally 
z = 3 in dilute solutions (non-draining) , but the case with 
z = (1 + 2v)/v (free draining or Rouse dynamics) is also 
conceptually important. The latter indeed corresponds 
to models in most (but not all) of numerical simulations. 
By denoting the translocation time r ~ Nq ^ our aim 
is to find these dynamical exponents a, /3 relevant to the 
driven translocation in terms of v and z. 

There is a growing evidence that a key element lies in 
the dynamical property of the polymer; it cannot respond 
to the driving force all at once. This becomes evident 
for rather weak forces / > /j ~ {k-oT / a)NQ^ , i.e., the 
thermal energy k^T divided by the equilibrium coil size 
i?o. In what follows, we make the length, force and time 
dimensionless by the units a, k^T/a and "qa^ /{k^T), re- 
spectively, with T] being the viscosity of the solvent. To 
analyze the case f > f^, let us see the sketches of a 
translocating polymer where the driving force acts on 
the segment at the pore {x = 0) as shown in Fig. 1. Note 
that X-axis is taken to be perpendicular to the wall with 
the pore. The one chain end is first sucked at time t = 0, 
while the polymer keeps the equilibrium conformation in 
X < 0. By the time t, M{t) segments have been already 
translocated, and the subsequent N{t) — M{t) segments 
in the anterior cis side (shaded in Fig. 1) are set in mo- 
tion under the tension, while the remaining Nq — N{t) 
segments in the rear part are still at rest. 



III. TWO PHASE PICTURE 

Although plausible, a mathematical formulation of the 
above picture is nontrivial. In fact, at any moment, each 
segment moves with its own velocity and the complete 
characterization of the transient process would require 
handling highly nonlinear equations in multi-dimensional 



where cr(x) is the line density of the segment and x = 
—R{t) is the location of the boundary (see Appendix II 
for this derivation). Note that we first consider only the 
tension-propagation stage (Fig. 1 (b), (c), (d)) since this 
is the dominant process giving the scaling exponents in 
the translocation. 

In addition, N{t)-th segment is just located in the 
boundary during the moving domain growth and then 
it still retains an initial rest position. Thus, by averag- 
ing over the initial conformation at t — 0, we find the 
relation 

N{ty^Rit). (2) 



IV. STEADY-STATE OF A DRAGGED 
POLYMER 

Here we write the dynamical equation of state relat- 
ing the steady-state velocity V, the extension R and the 
driving force / in the form; 



VR^Jif). 



(3) 



The scaling form J^{f) ~ and the segment line density 
at the boundary 



a{-R) ~ 

are expressed through the exponents: 
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As already mentioned, /j ~ Nq'^ signifies the onset of the 
nonequilibrium dynamics. Moreover, we find two char- 
acteristic forces /* ~ 1, /** ~ Nq. This leads to the 
identification of three distinct nonequilibrium regimes on 
the dynamical diagram as summarized in Fig. 2. The 
derivation and the physical meaning of eqs. ©-(jG]) and 
boundaries will be clarified in the following subsections 
concerning each regime. 
Trumpet (TP) regime: 

Let us first consider the situation with the driving force 
in the range /(i</</*~l. Here the blob concepts 
enable one to handle the space dependent length scale, 
whereby provide a generic description of the deformed 
dragged polymer irrespective of the model details [H, [l^l 
(Fig.[T](b)). Evaluating the dragging force and the elastic 
restoring force in terms of the blob size ^(x) ~ g{xY at 
each point, we have the local force balance equation as 
follows: 



(a) 



V{R 



(7) 



Substituting a; = combined with the total force balance 
/ ~ 1/^(0), we obtain the dynamical equation of state 
eq. ^ with eq. JS]) (top), which corresponds to a gen- 
eralization of the trumpet model for a tethered polymer 
under flow [l^l- At the boundary, there is the largest 
blob of size centered around x = —R + ^ij. Then, 
eq. ([7]) yields 



(8) 



which leads to the segment line density (7(- 
gn/^R =^ eq. dg (top) with eq. Q. 
Strong stretching (SS) regime: 

At the upper threshold of the Trumpet regime / = 
/* ~ 1, the number of segment in the smallest blob 
at the pore vicinity becomes very small ^(0) ~ 1 <J4> 
f(0) = ~ 1. For stronger forces, the blob descrip- 
tion is no longer valid, and one has to take account 
of the finite chain extensibility of the backbone (e.g.. 
Fig. [l](c), (d)). The line density such highly stretched 
parts is estimated as cr{x) ~ 1 at the scaling level, which 
indicates eq. ^ (bottom) with eq. Let us consider 
the limiting case of the strong driving so that the mov- 
ing domain of the chain is highly stretched close to the 
fuU extension R{t) ~ a{N{t) - M{t)) [21] in the entire 
stage of the translocation (Fig. [l](d)). This is realized 
when / > /** ~ Nq [27]. Neglecting the logarithmic 
correction due to the hydrodynamic interactions, the to- 
tal drag force is given by / ~ {N{t) — M{t))V , and then 
the scaling form of the dynamical equation of motion is 
specified by eq. ([5]) (bottom) irrespective of the solvent 
property (free-draining or non-draining). 

Stem-flower (SF) regime: 

In the intermediate between above two regimes /* < 
/ < /**, the moving domain takes a highly stretched 
conformation only in the early stage, for which the anal- 
ysis in the SS regime applies. This point will be discussed 
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FIG. 3: Plots of T{No,f) in the non-draining case (double 
logarithmic scale), (a) and the upper/lower plots in (b) cor- 
respond to the parameter sweepings indicated by arrows (i) 
and (ii)/(iii) in Fig. H 





Non-draining 


I' = 0.50 


f = 0.588 


EQ 


a = LOO, 13 = 1.00 


Q = 1.18, /? = 1.00 


TP 


a = 1.25, 13 = 0.50 


Q = 1.38, /? = 0.65 


SF 


a = 1.25, 13 = 0.50 


Q = 1.38, /? = 0.65 


SS 


Q = 1.50, 13 = 1.00 


Q = 1.59, /3 = 1.00 



TABLE I: Predicted exponents a and /3 in the non-draining 
case {z = 3). Note that 7 = v/a. 



in the section VI again. The moving domain grows with 
time, and the velocity V decreases towards ~ 1, at 
which the tension at the boundary becomes comparable 
to /*, that is a "Blooming point", before which the mov- 
ing domain composes of stem only. From then on, the 
conformation of the moving domain can be described as 
a stem-fiower (Fig. [D (c)). Compared to the SS regime, 
the chain in SF regime has smaller extension because 
of the coiling at the rear part. Still, scaling exponents 
associated with the entire domain in the SF regime is 
indistinguishable from that in the SS regime [2l| and 
satisfies Kit) ~ a[N{t) - M{t)] and / ~ V[N{t) - M{t)]. 
Consequently, we find eq. ([5]) (bottom); its distinction 
with that in the SS regime lies in a smaller numerical 
prefactor, while the conformation at the rear can be de- 
scribed by the TP model (eq. dH) (top)). Therefore, the 
SF regime shares (i) the dynamical equation of state with 
the SS regime, (ii) the segment line density at the bound- 
ary with the TP regime. 



V. DYNAMICAL SCALING: TENSION 
PROPAGATION 

A set of basic eqs. ([I])- (ID), combined with the appro- 
priate steady state conformational properties (eqs. ([5]) 
and ([HI)), leads to a differential equation: 

This should be solved with the initial condition at 
t = Tini, which plays a role in the correction to scaling (see 



Driving force / Number of segments 

FIG. 4: The same plots as those in Fig.[3]in the free-draining 
case. 





Free-draining 


V = 0.50 


V = 0.588 


EQ 


a = 1.50, /3 = 1.00 


a = 1.59, /3 = 1.00 


TP 


Q = 1.33, /3 = 1.33 


a = 1.43, /3 = 1.26 


SF 


Q = 1.33, fi = 0.66 


a = 1.43, /3 = 0.74 


SS 


a = 1.50, P = 1.00 


a = 1.59, /3 = 1.00 



TABLE II: The same as TABLE I in the free-draining case. 

below). Collecting the leading order terms, one obtains 
the dynamical scaling for the front evolution, which is 
asymptotically {t ^ Tini) valid, in the form 



/ = Tini 




FIG. 5: (Color Online) The sketch of the entire process in 
trumpet (TP) regime, (a) Initial equilibrium conformation, 
(b) Formation of the initial tensed blob, (c) Propagation 
stage of the tensile force retaining the TP shape, (d) / (e) 
Beginning/mid point of the post-propagation stage, (f) End 
of the translocation process. The tensed "moving domain" is 
shaded, while the translocated part is represented by a dotted 
curve. 



i?(i/)^i}, (10) 

where the growth exponent ^ = v j (\ ^ v — vq) charac- 
terizes the time evolution of the tension front and the 
force exponent /? — p(l ~ q) encodes the rescaling of time 
tj = tf^ with respect to the driving force. This stage per- 
sists until the tension reaches the rear end, from which 
the propagation time follows as 

Tp ^ Kf-^, (11) 

with a = v/j. Notably, such a scaling relation among the 
roughness, growth and dynamical exponents character- 
izing the nonequilibrium tension propagation process is 
reminiscent to that in the surface growth phenomena . 

Our main results are summarized in Figs. |3] and HI 
where we plot Tp{No, f) in the non-draining (z = 3) and 
the free-draining {z — {1 + 2v)/v) cases, respectively. 
These are expected to be valid for long chains, and in- 
deed capture a number of observations made in numeri- 
cal/real experiments. We also listed the concrete values 
of exponents for cases v ^ 1/2 and v — ^ 0.5876 
in tables 1 and 2. As a primal message, we stress that 
the identification of the appropriate regime (depending 
on iVo and /) as well as exponents z characterizing 
the equilibrium polymer coil is of a crucial importance 
for the proper comparison with experiments. For ex- 
ample, we first note that most simulations treat free- 
draining chains and the typical driving force is on the 
order of / ~ 1 (TP and/or SF regimes), and the re- 
ported values of a exponent are scattered around ~ 1.4 



in three dimension (y — 1^3). In particular, recent high- 
accuracy simulation results a — 1.41 ± 0.01 (Langevin 
dynamics) and a = 1.42 ± 0.01 (atomistic molecular dy- 
namic) [7] are in excellent agreement with our prediction 
a = 1.43. The change of a exponents with the increase 
in / and/or with the inclusion of hydrodynamic inter- 
actions has been examined in ref. Q, whose qualitative 
trend accords with ours. In experiments using double- 
stranded DNA 0, i|, the non-draini ng approximation 
with V = 1/2 would be appropriate [19[, and the typ- 
ical driving force range 10 pN would correspond to SF 
or even SS regime for short DNA chains. Then, the re- 
ported values a = 1.27 ± 0.03 {Sj and a ~ 1.4 [I are m 
good agreement with our prediction. On the other hand, 
much less attention has been paid to the /3 exponents. 
Usually, it has been implicitly assumed to be /3 = 1, but 
a recent study has claimed (3 ~ 0.8 for strong force |8j] in 
reasonable agreement with our prediction = 0.74 in SF 
regime {u = v^, free-draining). 



VI. FINITE SIZE EFFECTS 

How can we go one step further beyond the leading 
term scaling argument? Practically, resolving complica- 
tions associated with various finite-size and/or crossover 
effects is an important issue, for which we point out 
two main sources. The translocation time can be for- 
mally written as r = ri„i -I- Tp -I- Tpp. Figures 5, 6 and 
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FIG. 7: (Color Online) The sketch of the entire process in 
strong-stretching (SS) regime, (a) Initial equilibrium confor- 
mation, (b) Relaxation of first tensed segment, (c) Propaga- 
tion stage of the tensile force. (d)/(e) Beginning/mid point 
of the post-propagation stage, (f) End of the translocation 
process. 



FIG. 6: (Color Online) The sketch of the entire process in 
stem-flower (SF) regime, (a) Initial equilibrium conformation, 
(b) Relaxation of first tensed segment. (c)/(d) Propagation 
stage of the tensile force; (c) "blooming point" (end of the only 
stem conformation) / (d) stem-flower, (e) End of the propaga- 
tion stage (beginning of the post-propagation stage). (f)/(g) 
the post-propagation with (e) stem-flower/ (g) only stem, (h) 
End of the translocation process. The domains shaded by blue 
or orange correspond to flower or stem, respectively {x = —r 
is the junction point between stem and flower). 



with the exponent p given by eq. 

The point (ii) can be taken into account by the follow- 
ing estimates; 

f ?o - /-^ (/»</< /*) 

i?ini ^ I r(ri„i) ~ / (/*</< /**) (14) 

[ 1 (/** < /) 



7 show the full story of the driven translocation pro- 
cess in TP, SF and SS regimes, respectively. For long 
chains, the tension propagation stage (Fig. 5(b)— 7'(d), 
Fig. 6 (c)— 7-(e), Fig. 7(b)^'(d)) dominates the whole pro- 
cess in the nonequilibrium driven regimes (/ > /j), thus, 
Tp represents the translocation time r ~ Tp. For shorter 
chains, however, this might be severely affected by (i) 
the post-propagation stage (Tpp) and/or (ii) the effect 
associated with the initial conditions in eq. (|9]) (rini and 
i?ini). For the point (i), we can describe the process after 
t = Tini + Tp by setting dN/dt = in cq. ([T]) as follows: 

-R^^fP, (Tini + Tp < t < r) (12) 

where R{t) is the location of the rear end and ap- 
proaching to the hole (Fig. 5(d)^(f), Fig. 6(e)^-(h), 
Fig. 7(d)-^(f)) with Rin^i + Tp) Rq. The estimate 
for its characteristic duration time Tpp — t — (Tjni + Tp) 
is obtained by setting R{t) — 0, thus, 

Tpp Rlf-P ^ N^T" (13) 



f Co ^ f-' (/«</< /*) 

Tini ~ < G^i^S) f-PiSS) ^ Jl/^ (/*</< /**) (15) 

[ 1 (/** < /) 

The first one (/j < / < /*) comes from the formation 
of the initial tensed blob of size — (Fig- 5(b)), 
which becomes the microscopic time '~ 1 at the thresh- 
old f — f* ■ This microscopic time corresponds to Tini — 1 
in the SS regime as in Fig. 7(b). In the SF regime 
(/*</< f**)i we choose Tini to be the blooming 
time, i.e., the moment when the flower part starts to 
appear (Fig. 6(c)). Before that, there is only a stem 
part in the early period (Fig. 6 (b) — > (c)) during which 
the analysis for the SS regime applies after the micro- 
scopic relaxation time of first tensed segment t^^^^ ~ 1 as 
in Fig. 6 (b). These lead to the time Ti„i ~ G"(ss) /-/3(ss) 
after which the coiled flower- conformation appears at the 
rear part (a^ss), /^(ss) are the exponents for SS regime). 
Here, by denoting the stem length as r{t) (Fig. 6(d)), 
the pertinent segment number is G ~ [^'(Tini)]"'^^'^ from 
eq. where the stem length at the beginning of the 
late SF regime is r(Tini) — / (See Appendix III). Solving 
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eq. ([9]) with these initial conditions results in the shorten- 
ing of Tp. In particular, near the thresholds / — /**(— 0) 
and / — > /tt(+0), we find i?ini — > Ro, thus, Tp become 
vanishingly small. In the former, Tini is already close to 
Tp(/ — > /**(-l-0)). In the latter case, the contributions 
from Tini and Tpp are comparable, and Tpp crossovers to 
Tpp'^'' ~ Nq^~^^'^ , which is a result from the equilib- 
rium shape assumption [isj . These factors may result in 
the narrowing of the TP regime and, in general, cause 
the rounding-off the crossover at /j as well as /**. 

To summarize this section, one finds by comparing 
eqs. (|TT]) . (fT3)) and (fTSt that the propagation stage domi- 
nates, i.e., T ~ Tp for large N limit in driven translocation 
regime (/ > /j). For real systems with small N, however, 
the finite-size correction might be apparent. In partic- 
ular, Tpp may become dominant for short chains in two 
dimension, thus, an appreciable crossover scaling with re- 
spect to N would be expected as pointed out in ref. [l^ ■ 
On the other hand, for weaker forces (/ < /j), the equi- 
librium shape assumption would be valid, and t ~ ''pp'^^- 
We note, however, that in this weaker forces (/ < /j), the 
segment returning through the pore becomes apparent, 
and this may alter the exponent in Tpp*^^ , while such 
an effect is negligible under the strong driving force. 



VII. REMARKS AND PERSPECTIVE 

Once again, we emphasize that the key lies in the 
nonequilibrium transient response property inherent in 
a long fiexible chain. Our strategy to implement this 
is a two phase picture and applying a steady state ap- 
proximation for the moving domain. Before closing, let 
us make several remarks on the derived exponents, (i) 
Exponents in the translocation time a, /3 reflect the ex- 
ponent p, q characterizing the steady state conformation 
of the dragged chain (eqs. dS])-®). We have adopted 
eqs. (O, which is correct in the scaling limit, for 
the purpose of the general discussion, but a care should 
be taken in the comparison with real experiments. In 
particular, the logarithmic correction due to the hydro- 
dynamic interactions (in non-draining case) is neglected, 
and this might alter the effective exponent p, q , thus, 
a, (3, too, in the certain ran ge o f the interest j23l \2^ . 
Although the results in ref. (4 Il3l| are in good agreement 
with our prediction, further careful checks are awaited on 
it. (ii) As stated in Appendix I, the steady state approx- 
imation is by construction marginally valid to capture 
the salient features of the driven translocation, but there 
is some room to improve the theory by incorporating a 
slight deviation from the approximation, (iii) The pres- 
ence of several different regimes and the finite-size effects 
makes the accurate determination of the exponents very 
difficult (see also the remark (i), (ii) above). We expect, 
however, that it would be possible to capture a global 
trend in the exponents by carefully analyzing numeri- 
cal/experimental data obtained in widely different con- 



ditions, (iv) Aside from exponents, the proposed physi- 
cal picture is generic, and we believe that our formalism 
constitutes a minimal model, or a prototype in the prob- 
lem. From our point of view, it is more important to 
comprehend the underlying physics (even qualitatively), 
and seek for possible improvements of theory to describe 
the phenomena more adequately (see, for example, a re- 
cent attempt by Rowghanian and Grosberg [26]). (v) 
One of the interesting predictions concerns the beta ex- 
ponent; it is not necessarily unity in contrast to claims 
in most other works. We expect that such a trend (even 
qualitatively) could be detected by careful experiments; 
its reliable trend should serve as a stringent test for any 
theory of the driven translocation. 



VIII. SUMMARY 

In summary, we proposed dynamical scaling scenario 
applicable to the broad range of situations met in the 
study of polymer driven translocation. The paper was 
devoted to the detailed account of tension-propagation 
mechanism proposed in [3 its extension to cover 
the wider parameter space through the distinct nonequi- 
librium conformations. Also shown are the finite size ef- 
fect, various associated crossovers and their rounding off. 
We believe that these elements operate as a guide for 
further studies towards comprehensive understanding of 
the driven translocation. Last of all, we expect that the 
emergent physical picture may be relevant to many other 
polymeric systems driven by local and strong stimuli. 
Best examples would be found in the micro-manipulation 
experiments as well as the biopolymer functions in cells. 

Appendix I. — Validity of the steady-state ap- 
proximation: Let us look at the moving domain at 
some time during the tension propagating stage (e.g., the 
shaded part in Fig. 1 (b)). 

Since N{t)-th segment at the front have just taken part 
in the moving domain and the velocity profile would be 
smooth even around the front, there might be substan- 
tial velocity gradient there. If so, its magnitude can be 
estimated as 7^^, ~ V/^R, whose is the size of the 
largest blob at the front (eq. ([5])) and can be regarded as 
the boundary layer. Comparing the shear rate with the 
corresponding longest relaxation time t^^ ~ we find 
T^^7^^ — — 1 using eqs. Thus the concerned 

shear rate is marginal and one can apply the steady state 
approximation. 

As another possibility, one may suppose that an uni- 
form velocity gradient is entirely spread in the moving 
domain with the shear rate given by 'Jr ~ V/R. If so, 
the longest relaxation time in the entire tensed chain 
is tr ~ and we similarly have the relation as 

TR-yn ~ 1 [23 ■ Here tr was estimated as follows; The dy- 
namical equation of state ^ f — VRf / J{f) implies the 
spring constant fc/j = V f / J{f) and the viscous frictional 
coefficient Tr = Rf / J{f), which lead to the relaxation 
time tr ~ TR/kR = R/V. Although obtaining real ve- 
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locity profile in the moving domain requires a more thor- 
ough analysis, which is beyond the scope of the present 
study, one may expect it to be something between the 
above two extremes. We can therefore expect the steady 
state approximation would be valid to address the salient 
features in the driven translocation. 

Appendix II. — Mass conservation: Here we dis- 
cuss the derivation of eq. ([1]) from the integral form of 
the mass conservation in the moving domain; 



N{t) - M{t) = 



R{t) 



a{t, x)dx. 



(A-1) 



Taking account of the moving boundary R{t) associ- 
ated with the tension propagation, the time derivative 
of eq. (A-1) leads to 

mity^ . r ^d. + .(t,-i^)^ 

dt J-R(t) ot dt 

= j{t,-R) - j{t,{)) + a{t,-R)^, 

(A -2) 

where the continuity equation da(t^ x) / dt+dj {t, x)/dx = 
is used in the last equality. The pore current in 
eq. (A-2) would be strongly dependent on the pore speci- 
ficity, but can be eliminated by its definition j(t, 0) = 



dM/dt. Then, eq. ([T]) follows by assuming j{t, —R)(= 
a{t,-R)V{t,-R)) = a{t,-R)V{t), where V{t) is the 
steady state velocity of the moving domain. Under this 
assumption, the number of the translocated segments 
M{t) should be calculated from eq. (A-1), not from time 
integration of a{0)V{t). The latter quantity does not sat- 
isfy the conservation of mass as V{t) characterizing the 
chain deformation slightly differs from the velocity of the 
M{t)-th segment at the pore. We note that a better es- 
timate of the velocities V{t, —R) at the boundary (due 
to the slight derivation from the representative velocity 
V{t)) would improve the result. 

Appendix III. — Blooming point condition in 
the SF regime: At the end of the early SF regime (com- 
posed of stem only), the following conditions are met for 
the stem length R and the moving velocity 



riRV ~ /, f]aV 



(in the original unit). (A — 3) 



These are respectively the total force balance and the 
force balance at the rear end with the blob of size ~ a. 



From these, one obtains V 
dimensionless form V = V* 



1, R 



R ~ 

= 7" (Til 



the 



/■ 
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